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We investigate the effect of a homogeneous magnetic field on the thermal deconfinement transition 
of QCD in the large Nc limit. First we discuss how the critical temperature decreases due to the 
inclusion of Nf <C N c flavors of massless quarks. Then we study the equivalent correction in the 
presence of an external Abelian magnetic field. To leading order in Nf /N c , the fact that the decon- 
finement critical temperature decreases with the magnetic field depends solely on quarks behaving 
paramagnetically. Finally, we discuss the effects from a finite quark mass and its competition with 
magnetic effects. 



Introduction. The phase diagram of strong interactions 
in the presence of a classical, constant, and uniform mag- 
netic background has been attracting increasing interest 
in the last few years. Strong (Abelian) magnetic fields 
not only provide another control parameter to probe the 
phase structure of QCD but are also currently gener- 
ated in non-central ultrarelativistic heavy ion collisions at 
RHIC-BNL and at the LHC-CERN. In fact, these fields 
are believed to be the largest ever produced since the 
times of the electroweak phase transition in the early 
universe, reaching values on the order of B ~ 10 19 Gauss 
(eB 6 m 2 ) and even much higher Furthermore, lat- 
tice Monte Carlo simulations are not constrained by the 
sign problem in this case and can produce a trustworthy 
T x eB phase diagram, among other results. Neverthe- 
less, the mapping of this new phase diagram is still in 
its infancy and presents some conflicting pictures coming 
from different model calculations. 

In this Letter we study the behavior of the deconfining 
critical temperature T c in the presence of a strong mag- 
netic field in the large N c limit of QCD. This provides 
a well-defined setup for a clean, semi-quantitative de- 
scription by essentially counting powers of Nf/N c (with 
Nf being the number of quark flavors) when matching 
pressures for the confined and deconfined sectors. In 
the sequel we show that, in this limit, the deconfinement 
temperature decreases with the magnetic field for small 
Nf/N c in line with recent results from lattice QCD, pro- 
vided that quarks behave paramagnetically. We also dis- 
cuss how the critical temperature for the pure glue theory 
decreases due to the leading order correction in Nf/N c . 

All model calculations so far have suggested that suffi- 
ciently large magnetic fields, typically eB ~ 10 m 2 , could 
bring remarkable modifications in the QCD phase di- 
agram, from shifting the chiral and the deconfinement 
phase transition lines pl-fl^ to transforming the vacuum 
into a superconducting medium via p-meson condensa- 
tion [l3j ]. In particular, most model descriptions have 
predicted either an increase or a flat behavior for the de- 
confinement critical line as eB is increased to very large 



values. Exceptions can be found in Ref. [2J, where the 
critical temperature vanishes at a finite critical value of 
eB c ~ 25m 2 , featuring the disappearance of the confined 
phase at large magnetic fields, and in [3J, where vacuum 
corrections are disregarded, and T c diminishes with eB. 



The first pioneering lattice simulations 14], still with 
large values for the pion mass, also suggested a very mild 
increase of the critical temperature with eB. However, 
recent lattice simulations with physical masses [HI have 
shown that the critical temperature for deconfinement 
actually falls as the magnetic field increases. However, 
instead of falling with a rate that will bring it to zero at 
a given critical value of eB, it falls less and less rapidly, 
tending to saturate at large values of B in agreement with 
what one would expect from the phenomenon of magnetic 
catalysis 2, 17 1 . An exercise within the MIT bag model 
with the appropriate treatment of the subtleties of renor- 
malization at finite B has shown remarkable qualitative 
agreement with these lattice findings with respect to the 
behavior of T c (eB), i.e. it decreases and saturates for 
very large fields [18|. To the best of our knowledge, even 
if known to be crude in numerical precision and missing 
the correct nature of the (crossover) transition, this is the 
only description to date that captures the correct qualita- 
tive behavior of the deconfining transition in a magnetic 
background. 

Although a description of the deconfinement transition 
in the presence of an external magnetic field in terms of 
the MIT bag model is, of course, very simple, we believe 
it encodes the essential ingredient to provide a qualita- 
tive description of the behavior of T c x eB: confinement. 
The fact that the MIT bag model incorporates confine- 
ment (even if in its simplest fashion) seems to make it 
suited to describe the behavior of T c as a function of 
external parameters, as hinted by a previous successful 
description of the behavior of the critical temperature 
as a function of the pion mass and isospin chemical po- 
tential, as compared to lattice data, where chiral mod- 
els failed even qualitatively [lil, 0] • This suggests that 
confinement dynamics plays a central role in guiding the 
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functional behavior of T c and points towards a large N c 
description of the associated magnetic thermodynamics. 

Large N c thermodynamics. The large N c limit provides 
a great opportunity to study several aspects of QCD 
2l|-|2j| • Feynman diagrams are reorganized according to 
their dependence on N c and, when N c — > 00, only planar 
diagrams are relevant. The theory is still asymptotically 
free with a perturbative beta function defined in terms 
of the 't Hooft coupling A = g 2 N c and a renormalization 
group invariant energy scale Aqcd at which the associ- 
ated coupling becomes strong. While confinement has 
not been proven in this limit, it is widely believed that in 
the vacuum the physical degrees of freedom are weakly 
interacting (since interactions go as 1/N C ), colorless glue- 
balls. Nf quark degrees of freedom in the fundamental 
representation can be added to this theory and the cor- 
responding mesons are free when N c — > 00 while baryons 



see that there must be a deconfinement critical temper- 



become extremely heavy, Mb arvon ~ N c Aqcd 

Lattice QCD calculations |25[ show that the decon- 
finement phase transition of pure glue S U(Nr) gauge 
theory becomes first order when N c > 3 [26l429l | with 
a critical temperature lim jv c _> oo T c / ^fa — 0.5949(17) + 
0.458(18)/iV c 2 0, where a ~ (440 MeV) 2 is the string 
tension. The thermodynamic properties of pure glue do 
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not seem to change appreciably when N c > 3 
which suggests that large N c arguments may indeed cap- 
ture the main physical mechanism behind the deconfine- 
ment phase transition of QCD. 

The fact that limjv c ->oo T c/ '«/a ~ O(iV c ) and that the 
deconfining phase transition becomes strong first order 
can be readily understood using the following argument 
[HI . When N c — » 00 and the number of quark flavors Nf 
is finite such as Nf/N c <C 1, in the confined phase glue- 
balls and mesons are very weakly interacting and, since 
they are colorless, they only contribute to the pressure 
at 0(N®). Baryons are extremely heavy in this limit 
and do not contribute at zero baryon chemical poten- 
tial. Therefore, when N c — > 00 the only contribution to 
the pressure of the confined phase comes from the gluon 
condensate ~ AI^Aq CD , which we write in terms of the 
other renormalization group invariant energy scale a as 
Pconf — CgN^a 2 , where Co is a positive number of or- 
der 1. Moreover, it should be noticed that the entropy 
density in the confined phase vanishes. 

On the other hand, asymptotic freedom implies that 
in the planar limit the gluon pressure is P g i U on(T) — 
C SB fglue{T/yfo), where csb is a positive con- 
stant determined from the Stefan-Boltzmann limit and 
limr/^^oo f g iue(T/y/a) = 1. The function f g i ue de- 
pends implicitly on the 't Hooft coupling A(T) and, while 
its general form is not known when T ~ y/a, it should 
be a monotonically increasing function of T that interpo- 
lates from at T — » to 1 for T — > 00. Its form can be 
computed using perturbation theory at sufficiently high 
temperatures where A is small |33j . If Nf — 0, since the 
pressure is always continuous at any phase transition, we 



ature defined by the condition P g i ue {t<_ 
or, equivalently, 
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which implies that the solution Tj ^ is a pure number 
of 0(N°) that in general cannot be computed perturba- 
tively since it is obtained from the self-consistent equa- 
tion 
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Since f g i U e increases monotonically with T one obtains 



that T c (0) must increase with cq (note that the critical 
23l . |24| . temperature only vanishes if cq — > 0) (34| • Lattice calcu- 



lations have shown that Tc /^/a ~ 0.59 [30|. The phase 

transition to a Zn symmetric deconfined phase is then 

of first order when N c — > 00, Nf — 0, and the entropy 

density jumps from zero to a finite number of 0{N 2 ) at 
T (0)_ 

Leading Nf /N c corrections. The first correction to this 
picture appears with the inclusion of Nf flavors of mass- 
less quarks. The previous Zn c symmetry is broken ex- 
plicitly in the deconfined phase because of the presence 
of quarks. While the U(N f )®U(Nf) -> U(N f ) vector pat- 
tern of (spontaneous) symmetry breaking leads to N"j — 1 
Goldstone bosons (the "pions" ) , their contribution to the 
pressure of the confined phase is of 0(NjN®), being neg- 
ligible when Nf <C iV c . The quark contribution to the 
pressure of the deconfined phase is given by P qua rk (T) = 
N c N f T i c^ SB f quar k(T/y/a), where c qS B is the corre- 
sponding positive dimensionless number computed in the 
Stefan-Boltzmann limit and f qua rk is a monotonically in- 
creasing function of T such that lim T ; f qU ark {T) — 
1. 

When Nf/N c -C 1 the explicit breaking of Z^ c sym- 
metry is small, slightly smoothening the phase transi- 
tion into a very rapid crossover. The balance equation 
that defines the critical temperature T c modified by 
the quark flavors is obtained by equating the pressures 

Pconf = Pglue + P quark (T c (1) ) . Since f quark is 

a monotonic function of T, one should expect that the 
critical temperature gets shifted towards smaller values. 
In fact, in the limit where Nf/N c <C 1 one finds the self- 
consistent equation 
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It is possible to obtain the effect of the leading order 
Nf/N c correction on Tj 1 ' in terms of Tc°\ Keeping only 
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the first correction in Nf/N c , one may take Tc i-» 
inside the brackets in the equation above. Since the ratio 
f quark If glue is positive, one can define a new positive 
constant given by 



ci(N f ) = c 



1 Nf C q SB 



J quark I ^ J 



4iV c a 

C SB 



fgiue ( ) 
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Therefore, we see that the self-consistent equation for 
Tj 1 ^ has actually the same form as Eq. ^ and is given 
by 
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^ J glue 
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Thus, since c\(Nf) < cq and / 9 ; ue is monotonically in- 
creasing with T, we see that the leading effect of Nf 
massless flavors in the large N c limit is to decrease the 
critical temperature by a small amount of order Nf/N c 
with respect to TjP\ In other words, the addition of light 
quark flavors should decrease the value of the deconfinc- 
ment critical temperature, which has been observed in 
lattice simulations 35l- 37 1 . 

Large N c behavior of T c x (eB). The same line of ar- 
gument used above can be employed to show that the 
deconfinement critical temperature must decrease in the 
presence of an external magnetic field in the large N c 
limit of QCD. Assuming that Nf/N c <C 1 and m q = 0, 
the only contribution to the pressure of the confined 
phase that enters at 0(N 2 ) is given by the vacuum 
(B = 0) gluon condensate c%\N 2 o 2 . The gluon and 
quark condensates change in the presence of a mag- 
netic field [38|-|40j but these modifications are negligible 
in the large N c limit. The gluon contribution to the 
deconfined pressure is blind to the magnetic field, i.e., 

Pglue(T) = N^4 B fglue(T/V^)- 

On the other hand, the quark contribution P qua rk 
feels directly the effects of the magnetic field. In 
fact, the regularized contribution [41] of the mass- 
less quarks to the pressure is P, 

N C N, ■ 'KT.\rr4,„4 

flavors with electric charges {(N c — 1)/N C , — 1/N C } m 
units of the fundamental charge. Only the largest 
(~ N%) charge in each pair contributes to leading order 
in N f /N c . 

Notice that the function f qua rk is positive definite and 
must increase monotonically with T for a fixed value 
of eB until it goes to 1 in the high temperature limit 
T 3> y/<J, eB. Given our previous analysis for the case 
where Nf ^ and B = 0, one should expect that the 
critical temperature as a function of the magnetic field, 
T c (eB), must decrease with respect to by an amount 
of 0(N f /N c ). 
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2 ^ with 



pairs (Nf)/N c -C 1 being the number of pairs of quark 



This can be seen directly by equating the pressures at 



c N c a = N C T C c SB f g i ue [ -j= 



T eB 

+ N c N pairs (Nf ) T c 4 C q SB fquark (6) 



^T 2 e 

and noticing that, since the left-hand side of the equation 
above is fixed (and independent of B at this order), the 
addition of the quark contribution on the right-hand side 
must lead to a decrease of the critical temperature by 
an amount of order Nf/N c . In fact, the solution of the 
equation above for T c (eB), to leading order in Nf/N c , is 



T c {eB) 1/4 

. J glue 



where we defined 
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Since C2 (N pa i rs , eB) < cq, the same arguments used 
before guarantee that T c (eB) /T^ < 1 by an amount 
~ Nf/N c . Therefore, one concludes that, in the presence 
of an external magnetic field, the deconfinement critical 
temperature decreases with respect to its value for pure 
glue in the large N c limit of QCD. Whether T c (eB) is 
also lower than the critical temperature in the presence 
of Nf/N c flavors of massless quarks at B = 0, Tc , re- 
quires that C2 (Np a i rs , cB) < c\ . This implies that quarks 
should behave paramagnetically, with positive magneti- 
zation M (T c , eB) = dP quark (T Cl eB)/dB > for all val- 
ues of B. 

For a free gas implementation of the deconfined phase 
fgiue = 1 and, in the limit of strong magnetic fields 
eB/T 2 > 1, one finds that f quark ~ eB/T 2 0- Thus, in 
this case the magnetic suppression of the deconfinement 
critical temperature goes like eB N pa i rs / (N c a). In fact, 
this simple implementation in the limits of low and high 
magnetic fields provides a scenario in which the slope in 
T c (eB) decreases for large fields, as illustrated in Fig. 1. 

An eventual saturation of T r as a function of eB, as 



observed on the lattice [15( and in model calculations 
[lU , cannot be obtained using the limits discussed in this 
paper in a general fashion. The specific form of T c as a 
function of eB depends on the non-perturbative functions 
fgiue and fquark- In fact, in the large N c limit, our results 
indicate that T c (eB) can only be a flat curve if f qU ark is 
such that M(T c ,eB) is positive but it vanishes for large 
fields. In this scenario, the most likely explanation for 
the nearly flat curve found in the N c — Nf — 3 lattice 
study performed in Ref. [IB} is a net cancellation effect 
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FIG. 1: Cartoon of the T c x eB phase diagram in the large 
N c limit, using the approximation of free deconfined quarks 
and gluons. The numerical value 0.59 shown in the plot was 
extracted from Ref. 13011. 



that occurs for sufficiently large fields due to a magnetic 
field dependent contribution to the pressure below the 
phase transition. 

Quark mass effects. When m q ^ the pressure of the 
confined phase is increased by the quark contribution to 
the vacuum trace anomaly, N c Nfm q (—(qq)), where we 
used the fact that the quark condensate is negative. This 
is equivalent to a small positive shift of Cq, in leading 
order in Nf/N c : 



C 1-4 C„ 



= c 



1 Nf m q {—(qq)) 
4N C ^ 



(9) 



In the deconfined phase only the quark pressure will be 
affected by the quark mass effects, decreasing e.g. in 
perturbation theory [33|. In a large temperature ex- 
pansion, we may write: f quar k >-> fquark - c 3 m 2 /T 2 , 
where C3 is positive. Therefore, the critical temper- 
ature computation in this massive case follows the 
same steps that led to Eqs. (@| and (JSJ, with the 
substitutions cq i— > c mq > cq and f q 



f, 



quark 



cW q /T 2 



q > c 

^ fquark 



I — ^ 

m q — 

As a conse- 



quence, Ci{Nf,m q ) > ci{Nf,m q — 0) andTi'™"' 1 is higher 
than its massless counterpart, Tc. Interestingly enough, 
however, the corrections to Co and f quar k are respectively 
~ m q /^Ja and ~ {m q /Tc°^) 2 , being extremely small 
for reasonable values of quark masses, m q -C Tc . 
Therefore, in this large N c regime, we find that the criti- 
cal temperature as a function of m q is essentially flat, in 
line with lattice results [42l |43|. 

Of course, the explicit dependence of T c with respect to 
the quark mass (or equivalently the pion mass) will also 
depend on the details of the functions f g i ue and f quar k 



as well as the quark condensate. In the study performed 
111 

nana 

within an effective model implementation of 
the N c = 3 and Nf = 2 deconfined phase, T c /yfo was 
found to be nearly constant with respect to variations in 
the pion mass. 

In the presence of a magnetic field, the quark conden- 
sate and its influence on T c are unaltered at this order 
in Nf/N c , while the quark pressure receives magnetic 
contributions, becoming f q uark(T/\/a, m q /T, eB/T 2 ). 



Therefore, the critical temperature T ( 
tion of Eq. (JTJ) with C2 replaced by 
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While c mq > Co, which enhances the critical temperature, 
the brackets behave according to the quark pressure in 
the presence of B. Again, if quarks behave paramagnet- 
ically, fquark increases and there will be a competition 
between mass and magnetic effects, which may lead to 
a non-monotonic behavior of the critical temperature as 
the quark mass is varied. 

Final comments. It would be interesting to extend the 
discussion about the magnetic effects on the deconfinc- 
ment critical temperature to the Veneziano limit of QCD. 
In this case, one could also study whether chiral symme- 
try restoration coincides with the deconfinement transi- 
tion when Nf, N c — > 00 in the presence of an external 
magnetic field. 
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